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We perform an analytic and numerical study of parametric resonance in a three-neutrino frame- 
work for sub-GeV neutrinos which travel through a periodic density profile. Commensurate with 
the initial level of approximation, we develop a parametric resonance condition similar to the ex- 
act condition for two-neutrino systems. For a castle wall density profile, the v e — > oscillation 
probability is enhanced significantly and bounded by cos 2 #23- The CP phase 8 enters into the 
oscillation probability as a phase shift. For several cases, we examine the interplay between the 
characteristics of the castle wall profile and the CP phase and determine which profiles maximize 
the separation between oscillations with 8 — 0, ±?,7T. We also consider neutrinos which travel along 
a chord through the earth, passing from the mantle to core and back to mantle again. Significant 
enhancement of the oscillation probability is seen even in the case in which the neutrino energy is 
far from the MSW resonant energies. At 500 GeV, the difference between oscillation probabilities 
with 8 = and 8 = ^ is maximized. 
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INTRODUCTION 

The phenomenon of neutrino oscillations is a conse- 
quence of the fact neutrino weak interaction states are 
superpositions of the mass eigenstates. Relevant to os- 
cillation phenomenology are the neutrino mass-squared 
differences Ajk ■= — m 2 , and the PMNS mixing ma- 
trix which we denote as U [TJ [2] ■ The mixing matrix can 
be parametrized in terms of three real mixing angles 9jk, 
with j, k — 1,2,3 and j < k, and the Dirac CP phase 
6. The overwhelming majority of neutrino oscillation 
data fits quite well within this standard three-neutrino 
framework; however, when considering a subclass of the 
experiments, one can often accurately understand this re- 
stricted data in a two-neutrino framework, requiring only 
a single mass-squared difference Am 2 and mixing angle 
9. Solar neutrino experiments and long baseline (LBL) 
reactor experiments can be approximately parameterized 
by A21 and #12, and atmospheric and some accelerator 
neutrino experiments likewise can be described with A 32 
and #23- The ability to separate the data as such speaks 
to the smallness of the mixing angle #13 and the ratio of 
mass-squared differences |A 2 i/A 3 2|. 

Analogies exist between neutrino oscillations and me- 
chanical oscillations, and in particular, mechanical oscil- 
lators can exhibit large amplitude oscillations when some 
of the oscillation parameters change periodically. This 
parametric resonance is particularly prominent when the 
parameters change at twice the natural frequency of os- 
cillation. Examples of parametric resonance are pen- 
dula with vertically oscillating supports [3] and Faraday 
waves, surface instabilities created in a vertically oscil- 
lating container of fluid [3]. The possibility of paramet- 
rically enhanced neutrino oscillations was first noted in 
Refs. |6] ; it was shown that if neutrinos travel through 
matter with a particular periodic density profile the os- 



cillation probability can be considerably enhanced. 

As neutrinos travel through matter, the mixing an- 
gles and mass-squared differences are effectively modi- 
fied as described by the MSW effect |8] . Neutrinos 
which propagate long distances through matter of suf- 
ficient densities can incur significant interactions which 
are diagonal in flavor, as the interactions are mediated by 
the charged and neutral currents of the weak interaction. 
Neutral current interactions are democratic amongst the 
flavors, leaving the oscillation probabilities unchanged; 
however, since ordinary matter consists of electrons, pro- 
tons and neutrons, charged current interactions affect 
only the electron (anti-)neutrinos, modifying the oscilla- 
tion probability. The upshot is that the neutrino oscilla- 
tion parameters effectively change in a periodic manner 
if the neutrinos travel through matter with a periodic 
density profile, leading to possibility of parametric en- 
hancement. 

Parametric resonance in neutrino oscillations has been 
studied extensively through both analytical and numer- 
ical means [H |H1 I9H2"2"]. Given the small mass-squared 
differences of the neutrinos and the available energies 
from high flux sources, it is not possible to set up a 
tabletop experiment with the appropriate density pro- 
file so as to demonstrate parametric resonant oscillations; 
however, in Refs. [TDJ [TT], it was realized that the den- 
sity profile of earth's interior might provide a suitable 
laboratory. Indeed, the earth's density may be approx- 
imated as piecewise constant with two main regions-a 
mantle surrounding a denser core |23) . A periodic po- 
tential consisting of two piecewise constant regions of 
differing densities is often referred to as a "castle wall" 
potential. Exact analytic solutions for two-neutrino os- 
cillations through such castle wall profiles exist and serve 
as a fundamental tool for understanding parametric en- 
hancement for core-crossing trajectories [51[T31[T5]. Spe- 
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cific applications consider atmospheric neutrinos which 
travel through the earth (T4j PT8H22] . Relatively exhaus- 
tive semi-analytic and numerical studies for neutrino os- 
cillations in the earth were done in Refs. [JJ|] [20] where 
resonance regions are shown to follow from generalized 
amplitude and phase conditions. For semi-analytic treat- 
ments of three-neutrino oscillations, one typically reduces 
the problem to an effective two-neutrino system at vary- 
ing levels of approximation. In this manner, one may 
incorporate the Dirac CP phase into the analysis, some- 
thing not possible in a pure two-neutrino theory. This is 
what is done in Refs. [201 HI] where the consequences of 
CP violation to the oscillation probability is considered 
for neutrinos with energies in excess of 1 GeV traveling 
through the earth; particular attention is paid to the in- 
terference between oscillations due to the A 2 i and A 3 i 
mass-squared differences. 

We will examine, herein, the interplay between CP 
violating effects and parametric resonance for sub-GeV 
neutrinos. Using approximations relevant for sub-GeV 
neutrinos and mantle/core densities, the three neutrino 
system can be cast into an effective two-neutrino system 
via the so-called propagation basis [53]. We apply ex- 
isting work on two-neutrino parametric resonance to a 
novel semi-analytic study of sub-GeV neutrinos, includ- 
ing CP violation. At a level of approximation commensu- 
rate with that used to rotate to the propagation basis, we 
find a condition for parametric resonance similar to the 
two neutrino case, and upon implementing this condition, 
we are able to assess the value of the u e — > oscillation 
probability at the end of an integer number of periods of 
the matter potential. We show that V efl is enhanced and 
bounded by cos 2 #23 here. Also, we are able to determine 
the characteristics of a castle wall profile that will lead 
to maximum separation between "P e/J curves employing 
different values of the CP phase <5 = 0, ±f,7T. We then 
turn briefly to the situation in which neutrinos do not 
travel through an integer number of periods as this is 
relevant for atmospheric neutrinos passing through the 
earth. Again, we focus upon parameters which imple- 
ment parametric resonance and then look at oscillations 
with different values of the CP phase. 

Our focus upon sub-GeV neutrinos is motivated by 
previous work which has shown that for long baselines 
the v e — > oscillations driven by the solar mass-squared 
difference contain relatively sizable terms proportional to 
sin which are CP odd [24 27 . As we enter an era 
of precision neutrino experiments, evidence for a nonzero 
value of mounts and with it the possibility of measur- 
ing the level of CP violation, if any, in the neutrino sector. 
Strict upper bounds on the magnitude of #13 were ini- 
tially established by the CHOOZ reactor experiment [28] : 
however, a recent reevaluation of the reactor neutrino 
flux [21] has somewhat relaxed this upper bound. Fur- 
thermore, hints of nonzero #13 come from joint analyses 
of solar neutrino and KamLAND data [501 - I52"] . Though 



statistically less significant, analyses of atmospheric neu- 
trino experiments also favor a nonzero reactor mixing 
angle |33t435] . Accelerator — ¥ v e appearance exper- 
iments MINOS [3B] and T2K [37] have both detected 
electron neutrinos above the expected background, fur- 
ther evidence for non-zero 6*i 3 . A global analysis of this 
neutrino data, excluding recent reactor experiments, in- 
dicates a value of #13 differing from zero by more than 3<r 
[55] . Perhaps most significant are the data from two re- 
actor D e disappearance experiments; both Daya Bay |39j 
and RENO [30J report nonzero values of sin 2 26*13 at the 
5-cr level. 



OSCILLATION IN MATTER 

The ultrarelativistic limit of the evolution equation for 
a neutrino of energy E is 



(1) 



where we define the matrix M. = diag(0, A21, A31). We 
employ the parametrization used in Ref. [24j 



U = U 1 (023)D S U 2 (O 13 )U 3 (O 12 ), 



(2) 



where Uj(9) is a proper rotation by angle 8 about the 
j-th axis and D$ = diag(l, l,e lS ); this is different from, 
but equivalent to, the standard parametrization found in 
Ref. [31]. 

When neutrinos travel through matter, the Hamilto- 
nian accrues an effective potential due to the coherent 
forward scattering of the neutrinos upon electrons, pro- 
tons, and neutrons which comprise the matter [7J [5] • We 
include this effective potential in the evolution equation 



id t v 



2E 



UMlP +V{x) 



(3) 



Neglecting the (irrelevant for oscillations) neutral cur- 
rent interaction, the operator V(x) exclusively acts on 
the electron flavor with a mag nitude V = V2G F N e (x), 
where Gf is the Fermi coupling constant and N e is the 
electron number density. We note that for anti-neutrinos, 
we need to change the algebraic sign of this potential and 
the CP phase 5. We shall consider only neutrinos below. 

For sub-GeV neutrinos traversing the earth, matter ef- 
fects are most easily addressed in the propagation basis 
developed in Ref. [24] . We will briefly review this deriva- 
tion. As the Ui(023) portion of the mixing matrix com- 
mutes with V, we may rewrite the evolution equation 



1 

2E 



U 3 {9 l2 )MU 3 {9 12 ) ] + U 2 (9 13 )^VU 2 (9 13 ) 



with v' = U 2 (9 13 )WtU 1 (9 23 )^. 



v 

(4) 
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By conjugating the Hamiltonian in this basis via a lo- 
cally defined U 2 (6), this new propagation basis can be 
approximately described by a Hamiltonian H which is 
block diagonal. This correction to #13 is given by 



tan 29 



2sin26 lz EV 



A31 - s? 2 A 2 i - 2cos20 l3 EV 



(5) 



where we use the shorthand s\2 := sin #12. The density of 
the earth's interior has an upper bound around 13 g/cm 3 
|23) : this results in a maximum effective potential V ~ 
5 x 1CT 13 eV. As A31 ~ 2.4 x 1CT 3 eV 2 , the mass-squared 
difference is the dominant term in the denominator of 
Eq. (|]); for E ~ 1 GeV, one has e := 2#V7A 3 i < 0.4. 
For sub- GeV energies, one may approximate Eq. ^ as 



sin26»i 3 £;^ 



A 



31 



This correction results in a modified mixing angle 



(6) 



(7) 



With this additional rotation, we define locally the prop- 
agation basis with v — XJiiff^v' and Hamiltonian H 



H 



H 

A 31 /2£- 



where the block is given by 



H 



1 

2E 



'12 



A 2 1 

Sl2 



2EV S12C12A21 
C12A21 c 2 2 A 2 i 



13 



(8) 



(9) 



Through the definition of 1/ , we directly relate the prop- 
agation basis to the flavor basis via v = tftv where 



U 



u 1 (e 23 )D s u 2 ( 



Cl3 
-Sl 3 S23e 
-S13C23C 1 



iS 





C23 



Sl3 
Cl3S23e 



-S23 Ci 3 C 2 3e 



id 



(10) 

(11) 



Thus, electron and muon neutrinos can be written in the 
local propagation basis as 




"Sl3S23e 

C23 



Cl3S23e 



-t<5 



(12) 

As the correction 6 depends upon the local density, 
we must consider its temporal (spatial) derivative in the 
evolution equation 



d t v = t [W)V + U 2 (9)%i/. 



(13) 



Letting A2 be the generator of the rotation so that 
U 2 {6) = e iex \ we have 



d t U 2 {ey = -i\ 2 U 2 {9)^d t t 



(14) 



Dropping insignificant terms, one may write the evolu- 
tion equation in the propagation basis as 



id t v = {H + \ 2 d t 6)v. 



(15) 



Considering only propagation through matter of con- 
stant density, the term dtO vanishes, and our evolution 
equation is 



idtv = HP. 



(16) 



The block H in this Hamiltonian can be easily diago- 
nalized in closed form with eigenvalues X±. Of dynami- 
cal relevance is the difference in these eigenvalues which 
yields the effective constant density mass-squared differ- 



A^ = A 2 i\/cos 2 26*12(1 - E/E R ) 2 + sin 26 12 , (17) 



where we have defined the resonance energy to be 



Er — 



A 2 1 COS 2#i2 

2Vc 2 3 ' 



(18) 



Fixing the solar mixing angle # 12 = 0.58, we find the 
resonance energy in the mantle of density 4.5 g/cm 3 to be 
En ~ 100 MeV; in the core of density p = 11.5 g/cm 3 , the 
value is En ~ 40 MeV. The mixing angle which achieves 
this diagonalization satisfies 



sin 26>™ 



sin 26-t 



cos 2 26>i2 (1 - E/E R ) 2 + sin 2 26 x 



(19) 



At resonant energy, the effective mixing angle in mat- 
ter, 6 |m ! , results in maximal mixing for these two neu- 
trino states in the propagation basis; this is termed the 
MSW resonance. Additionally, matter effects require 
an accommodation to the other mass-squared difference 
= A 3 i — 2i?A_, though this correction is dominated 
by the vacuum value of the mass-squared difference. 

In the analytic work that follows, we will be primarily 
interested in the oscillatory region for the small mass- 
squared difference A 2 i. We will assume that the oscil- 
lations due to the two larger mass-squared differences 
cannot be resolved at the baselines of interest L; that is, 
we will take 



sm 



A31L 
AE 
A31L 
AE 



A 32 £ 
AE 

A 32 L 
AE 



_ 1 
~ 2 

= 0. 



(20) 
(21) 



The upshot is that for sub- GeV neutrinos traveling 
through the earth the propagation basis provides us with 
a density-dependent effective two neutrino framework. 
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TWO FLAVOR PARAMETRIC RESONANCE 

We will review parametric resonance within the con- 
text of a pure two neutrino system, say, v e and v^, and 
briefly rehash known results. This construction can be 
suitably adapted to describe parametric resonance in an 
effective two-neutrino framework for sub-GeV neutrinos 
traveling though matter of terrestrial densities. We will 
denote the lone mixing angle as 9 and mass-squared dif- 
ference A. An exact solution for two neutrinos traveling 
through a castle wall potential was developed in Ref. [6 
and expounded upon in Ref. [T3"] . 

Following Ref. [TS], we will review the exact solution 
for neutrinos traversing a general periodic potential and 
then specify to the castle wall solution. Without loss of 
generality, one may choose the neutrino Hamiltonian to 
be traceless. If the two-neutrino Hamiltonian H is not 
traceless from the start, one may add to the Hamiltonian 
with impunity any multiple of the identity, in particu- 
lar — |tr(JJ)I; upon solving for the time evolution of the 
system, such multiples result in an immeasurable overall 
phase. Thus, we may take the Hamiltonian to be of the 
form 



H = 



—a(x) fi(x) 
(3(x) a(x) 



(22) 



where the real functions a(x) and (3(x) may depend on 
position by virtue of their density dependence. If the 
Hamiltonian is expressed in the flavor basis, then these 
functions are 

a(x) = ~cos(26)-~V(x), p = ~sm(29). (23) 

Denoting the period of the Hamiltonian as L, we have 
H(x) = H(x + L). 

As the Hamiltonian is Hcrmitian, the evolution of the 
system is unitary v(x) = U(x)v(Q). Given this, we may 
use the Pauli matrices to write evolution through one 
period 

U(L) = Y - ia ■ X; (24) 
unitarity demands that the real quantities satisfy 



Y 



1. 



This can be written in terms of a phase 4> 
U{L) = exp \-i(a ■ X)$ 

with unit vector X = X/|X| and 

cos$ = Y, sin$ = |X| 



(25) 



(26) 



(27) 



In this last formulation, it is quite easy to see that, after 
k periods, the evolution operator can be written as 



U(kL) = [U(L)] k = exp -i(a ■ X)fc$ 



(28) 



Thus, if the neutrino state is initially i/(0) = v e 
then after k periods the state of the system is 



v{kL) = 



cos(/c<J>) - iX 3 sin(fc<&) 
(X 2 - iX x ) sin(fc$) 



(29) 



As X is a unit vector by definition, we can redefine the 
terms which involve X\ and X2 by introducing a phase 
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X x + i X 2 = e 17 \l 1 - XI , (30) 



so that we may rewrite the neutrino at baseline kL as 
cos(A;$) - iX 3 sin(/c<I>) \ 



v(kL) 



-Hi 



*7 



1 - Xf sin(fc$) 



(31) 



In this form, it is clear that a maximum oscillation v e — > 
Vn can be achieved if X 3 = 0. This is the condition for 
parametric resonance for a general periodic Hamiltonian. 

An analytical expression for A3 is hard to come by for 
a general density profile; however, a tractable solution 
does exist for the castle wall potential Q21 [15] . Explic- 
itly, the castle wall potential is defined as the periodic 
piecewise-constant function given by 



V(x) 



Va for Q < x < L A 
Vb for La < x < Li 



(32) 



with the periodicity condition V(x + L) = V(x) where 
L = La + Lb ■ 

Within one of the constant density regions, the effec- 
tive mass-squared difference ^ and mixing angle 9a,b 
in matter can be determined by diagonalizing the Hamil- 
tonian as in Eqs. ( 17|19 1. The evolution operator through 
one period is composed of the constant density evolution 
operators U A,B (x) 



U{L)=U B {L B )U A {L A ) 



(33) 



where the constant density operators can be expressed as 



U A {L A ) 



c A + ic 2 e A SA -is26 A SA 
-is 2 e A SA c A - ic 2 e A SA 



(34) 



with the dynamic terms defined to be ca = cos tpA and 
sa = sint^ where ipA = AaLa/4:E. An analogous ex- 
pression exists for 14 (Lb)- Using the properties of Pauli 



matrices, one can express U(L) in the form of Eq. (24) 
0Q3] with 



Y = c A c B - (s 2 e A s 2 g B + c 2 o a c 2 o b )s a sb (35) 

/ S A CBS 2 A + S B CAS 2 9 B \ 
X = I (s 2 g B C 2 A - S 2 g A C 2 g B )sASB \ ■ (36) 
-S A C B C 2 g A - 8 B CAC 2 g B 



Thus, the condition for parametric resonance in a castle 
wall potential is 



SACBC 2 g A 



SBCAC 2 g E 



0. 



(37) 
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Baseline (L) 

FIG. 1: Oscillation probability — > through a castle wall 
potential using the following input: 9 — 0.1, A = 7.6 x 10~ J 
eV 2 , E = 200 MeV, p A = g/cm 3 , L A = 3255 km, p B = 10 
g/cm 3 , L B = 3212 km. 




Baseline (L) 

FIG. 2: Oscillation probability ^ e — s> ^ M through a castle wall 
potential using the following input: 9 — 0.1, A = 7.6 x 10 -5 
eV 2 , E = 200 MeV, p A = g/cm 3 , L A = 4882 km, p B = 10 
g/cm 3 , L s = 4819 km. 



Parametric resonance can be achieved via the "half- 
wavelength condition" in which L A and Lb are equal to 
an integer plus one-half (local) oscillation wavelengths; 
this amounts to c A = = c B - Such a half- wavelength 
scenario is pictured in Fig. [I] with a vacuum mixing an- 
gle of 9 = 0.1. With this mixing angle, the maximum 
vacuum oscillation probability for appearance would 
be 0.04, yet after four periods of the castle wall potential, 
the probability is unity. In general, the oscillation does 
not attain unity at the end of a period but, rather, at 
some point in between. For a profile satisfying the half- 
wavelength condition, one may determine from the defi- 



nition of Eq. (27), and the expression for Y, Eq. (35), 



that the oscillation probability will be unity at the end 
of the fcth period in the event that there exists an integer 
n such that 



(38) 



2k\e A -e B \ = [n+- i tt. 



For the parameters used to generate Fig. [T] it just so 
happens that the above is approximately satisfied for k — 
4, i.e., 8\9 A - B \ = 3.49 tt. 

More generally, the condition for parametric resonance 
can be satisfied whenever 



ta.mp B 



C 29 A 
C 26 B 



tan Lp a- 



(39) 



In Fig. [2] we demonstrate such a scenario. Of note in 
this example is the fact that the oscillation probability 
reaches unity roughly halfway between the fifth and sixth 
period of the matter potential. 



THREE FLAVOR PARAMETRIC RESONANCE 

Considering all three flavors, we can now study 
parametric resonance for sub-GeV neutrinos traveling 
through a castle wall potential for densities less than 15 
g/cm 3 . We saw above that for sub-GeV neutrinos the rel- 
evant oscillations can be cast, to an good approximation, 
in the form of two-neutrino oscillations after rotating to 
the propagation basis. As the rotation to the propaga- 
tion basis is density dependent, we must use a transi- 
tion matrix at the boundaries of the regions of constant 
density. In our analytic treatment, we will make use of 
this and further approximations; however, our numeri- 
cal computations will model a full three-neutrino system 
sans approximation. 

We consider the same castle wall potential as in 
Eq. (32). Using the sub-GeV approximation, the neu- 



trino state after one period is given by 



v(L) = U B e 



— iHs Lr 



U 2 {9 1 



u\m (40) 



where U a,b is evaluated for the constant potential Vab 
and 9a,b represent the matter corrections to #13, Eq. (|5j), 
for Va,b- To leading order in £#13, the transition matrix 
between the two regions is 



u 2 {e B -e / 



29 13 E(V B - V A ) 
A31 



A, 



(41) 



The evolution operator in the propagation basis within a 
constant density region is 



U ee (x) U ell {x) 



U A {x) := e-^ AX = \ U, e {x) U^{x) | ; (42) 

Urr(x) 



() 



we make an analogous definition for U B (x) = e lffflI . 
After one period, the neutrino state is to leading order 



v{L) = U B U{L)U\v{0) + 0(6 13 e) 



(43) 



where the leading order contribution to the evolution op- 
erator is denoted by U(L) = U B (Lb)U a (La), consistent 
with Eq. ( 33 1 . This is the dominant contribution to sub- 



GeV oscillations in the earth, but we examine the 0(e# 13 ) 
correction. Returning to the transition matrix, Eq. (41 1, 
we note 



BjiA 

ee w rr 

/-.', r ,x ;-/ 1, i i—l n () U B e U A 

—U B T U A e —U B T U A ^ i ) 



o u u, L ;u 

U"(L b )X 2 U a {L a ) = I 



(44) 

where U^p = U^^La), etc. Supposing v(0) = v e , the 
probability of detecting is 



V m {L) = \si 3B S23Cl3 A e lS U ee (L) - C 2 3Ci3 A U fie (L)\ 2 

,-2 ~2 2 
+ s 13a C 13 b s 23 

+(6 B - e A ) 2 c 2 S 2 3 d 2 \U ee (L)\ 2 , (45) 



assuming the "atmospheric" oscillations average to zero 
(U TT ) ~ and neglecting terms 0(df 3 ). The term in 
the oscillation probability proportional to (Ob — 6a) 2 is 
the leading order remnant from the transition matrix be- 
tween the constant density regions, and its size is order 
0(6\ 3 € 2 ). Neglecting this term will result in a disconti- 
nuity in the oscillation probability across the boundary 
between regions; however, the 0(e 2 ) factor suppresses 
the significance of this discontinuity. We shall neglect 
this term in our semi-analytic analysis; this amounts to 
setting the transition matrix to the identity. 



With this simplification, the neutrino system in the 
propagation basis effectively consists of two states as the 
evolution of the v r state decouples. The analysis of two- 
neutrino parametric resonance can be carried over whole- 
sale with one adjustment for the effective potential in the 
propagation basis, namely, Va,b ^ c 2 3a b Va,b- Given 
this, if one begins with an initial state i/(0) = v e , then 
the neutrino state in the propagation basis at x = kL is 



v[kL) 



( ci3 A [cos(fc$) - iX 3 sin(fc$)] \ 
di3 A H^Jl-X 2 S m(k^>)} 



\ 



(46) 



I 



Si3 A l^T T {kL) 

The probability for a detection after k periods is 



V eii (kL) 



-13, 



si3s s 23[cos(fc$) — iX 3 sin(fc$)] 

2 



1 - X 2 sin(fc$) 



+ c 13 f? S 



13a ,5 23- 



(47) 



This expression for the oscillation probability is con- 
siderably more complicated than its purely two-neutrino 
analogue. In particular, for arbitrary vacuum mixing an- 
gles and mass-squared differences, it is clear that the ap- 
pearance oscillation probability v e — > cannot generally 
become unity via parametric resonance. This is no sur- 
prise, given the additional v T oscillation channel. Still, 
the question remains as to how a particular matter profile 
might maximally enhance the oscillation through para- 
metric resonance. To develop a parametric resonance 
condition, it is best to examine the oscillation probabil- 
ity to leading order in 9\3 



V^(kL) =sin 2 (fc$) C 2 3 (l-X 3 2 ) 



1 - 



2#i 3 S 2 3 cos(7 - 5)X 3 



C23 



i-xl 



2 sin(fc$) cos(fc$)0i 3 s 23 c 2 3 sin( 7 - S) yf 1 - Af + 0(8{ 3 ) 



(48) 



Overall, all terms in Eq. (48) are modulated by (at least) 
one factor of \ 1 



A|, so that the condition for para- 
metric resonance will be a perturbation of the purely two- 
neutrino condition, X 3 = 0. In fact, one can show that 
the value of A 3 which maximizes V e u is of order O(0\ 3 ). 



The term in Eq. ( 48 ) which is linear in X3 also has an 
explicit factor of #13. As our approximation is only valid 
up to 0(613), we must, for consistency's sake, effectively 
adopt the purely two-neutrino condition for parametric 
resonance, A3 = 0. Combining the remaining terms, we 
find, consistent to our level of approximation, the oscil- 



lation probability to be 

V e)1 {kL) = c 2 3 sin 2 [/c$ + ip] 



(49) 



with ip = —6*i3 tan 6*23 sin(7 — i5) when A3 = 0. Near the 
boundary of an integer number of periods, we see that 
the oscillation probability is bounded by c 2 3 rather than 
the unit bound in the purely two-neutrino framework. 
Also, terms linear in #13 which are CP odd, enter only as 
a phase shift in the oscillation probability. 

In this three neutrino system, the bound of C23 in 



Eq. (49) arises from the projection to the propagation 



basis which permitted the use of the two-neutrino anal- 
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ysis. If we work within the effective two-neutrino pic- 
ture in the propagation basis, then the amplitude of the 
oscillation probability is sin 2 29\2 ; as with the two neu- 
trino case, parametric resonance allows one to saturate 
the transition probability at unity. Returning to the fla- 
vor basis, we then expect the probability to saturate at 



C23 . It is useful to compare Eq. ( 49 ) with the oscillation 
probability for sub-GeV neutrinos in matter of constant 
density. Referring to Ref. [57], let us only consider the 
v e v u oscillation probability for the situation in which 
CP is maximally violated with 5 = |; for other values 
of (5, similar arguments hold. To leading order in #13, we 
have 



sin 2 2#i2C23 sin 2 



AE 



(50) 



with the phase <fi ~ — 0i3tan023/sin2#i2. Via para- 
metric resonance we can saturate the sin 2 2#i2 bound; 
sending sin 2 26* 12 H> 1, we find an expression similar to 



Eq. (49) 



solid [black] curve the effective two-neutrino approxima- 
tion in which we take as the identity the transition matrix 
between boundary layers. This curve is superimposed 
upon the results of an exact three-neutrino computation, 
plotted as the dashed [red] curve. For the three-neutrino 
curve, we average over the A31 oscillations so as to mimic 
a detector's finite energy resolution; remnants of these os- 
cillations appear as the higher frequency wiggles in the 
curve. The effective two-neutrino approximation accu- 
rately captures the oscillations driven by the A 2 i mass- 
squared difference. There is a slight discontinuity in the 
solid curve at the boundary between regions A and B; 
however, it is not too severe as the product C$13 is rather 
small. Finally, keeping only terms linear in #13, we were 
able to determine the oscillation probability after an inte- 



ger number of periods, Eq. (49 ). In the figure, we plot the 
value of Eq. (49) using the cross, X. This agrees rather 
well with the other two curves though it does have a sys- 
tematically lower value than the more exact treatments; 
we trace this to a positive term of order O(0\ 3 ) that has 
been neglected. 




2L 3L 4L 

Baseline (L) 



5L 



6L 



FIG. 3: (Color online) Oscillation probability v a — > 
through a castle wall potential using the following parameters: 
0i2 = 0.58, 0i3 = 0.15, 23 = 0.785, 5 = 0, A 2 i = 7.6 x 10~ 5 
eV 2 , A 3 i = 2.4 x 10~ 3 eV 2 , E = 200 MeV, p A = 4.5 g/cm 3 , 
L A = 3161 km, p B = 11.5 g/cm 3 , L B = 1597 km. The dashed 
[red] curve uses the exact three-neutrino framework, averag- 
ing over the A31 oscillations. The solid [black] curve employs 
the effective two-neutrino approximation in which the tran- 
sition matrix between regions A and B is taken to be the 
identity. The cross, X , represents the value of Eq. ( 49 1 at 
points kL. 

In Fig. [3j we compare the approximate analytic treat- 
ment for neutrinos traveling through a castle wall po- 
tential with exact numerical results. Realistic values of 
the mixing angles, mass-squared differences, and den- 
sities have been chosen so as to satisfy the half wave- 
length condition. The parameters which are germane 
to our approximations have the values #13 = 0.15 and 
e := 2£ , V B /A 3 i = 0.07. In the figure, we plot as the 




2L 3L 
Baseline (L) 

FIG. 4: (Color online) Oscillation probability v e — > 
through a castle wall potential using the following common 
input: 0i2 = 0.58, i3 = 0.3 , 23 = 0.785, A 2i = 7.6 x 10" 5 
eV 2 , A31 = 2.4 x 10~ 3 eV 2 , E = 200 MeV, p A = 4.5 g/cm 3 , 
L A = 3239 km, p B = 11.5 g/cm 3 , L B = 1707 km. The solid 
[black] curve has 5 = 0; the dashed [red] curve has 5 — n; the 
dotted [green] curve has 5 = ? ; and the + [blue] curve has 
5 = ^? . The A31 oscillations are averaged over. 



In Fig. [4J we implement parametric resonance, X3 = 0, 
via the half-wavelength condition, a three-neutrino ana- 
log to Fig. [TJ For the neutrino oscillation parameters, we 
use the best fit values from the global analysis of oscil- 
lation data in Ref. [35], save #13. In order to accentuate 
the effects of this mixing angle, we set #13 = 0.3 which 
is roughly 10-cr larger than its best fit value [38]. For 
this matter profile, the oscillation phase for a single pe- 
riod is $ = 2.49, and the half-wavelength condition forces 
X2 = 1 so that 7 = f ■ With maximal mixing for #23, 
the phase offset for the oscillations is ip = — 0.3cos<5. De- 
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spite the fact that this is relatively small, \ip\ < 0.3, the 
phase can have a large impact as to where the (absolute) 
maximum oscillation probability occurs, since Eq. ( |49[ ) is 
a function of a discrete number of periods kL. For four 
different values of 8, we see widely varied traces for V eit . 
The solid [black] curve plots the oscillation probability 
for the CP conserving case of 8 = 0. Focusing only upon 
the curves at points kL, this curve attains is maximum 
value near x « 2L, as 2$ + ijj ss As a contrast, the 

dashed [red] curve also plots CP conserving oscillations, 
but with 8 = ir. Relative to the first case, the sign of 
-0 changes, and we find the maximum value of V efJ , is at- 
tained for isi3L given that 3& + ip ~ 2.5-7T. We also plot 
two cases of maximal CP violation. The dotted [green] 



curve has 5 = ^, and the + [blue] curve has 8 — 



3tt 
2 



As tp = for both of these cases, they both intersect for 
integer multiples of the period. Also, at these points, the 
oscillation probability takes a value intermediate of the 
two CP conserving cases. 

More notable, perhaps, is the separation between the 
CP violating and CP conserving curves at the end of 
each period. For a general matter profile, the difference 
between the CP conserving curves, i.e., 8 = and 8 = n 
[or equivalently ±#13] at the point kL is 



-#13 sin 2023 sin 7 sin 2fc$. (51) 



The maximum separation between the two curves occurs 



when fc$ = 



f for some integer n and 7 



7T 3_7T 
2 ' 2 



|7Wo-7WJ < 1 013 sin 29 



23 



(52) 



These conditions can be trivially satisfied with the half- 
wavelength condition. Incidentally, for sub-GeV neutri- 
nos propagating through a constant density region in the 
earth's mantle or core, it was shown in Ref. 142 that 



< |0i3 sin 2023 sin 40 



23 > 



(53) 



For a mantle density of p = 4.5 g/cm 3 , the difference 
between the oscillation probabilities in constant density 
matter for ±0i3 is suppressed by a factor of | sin 40™2 1 = 
0.75 relative to the castle wall profile, whereas for a den- 
sity of 11.5 g/cm 3 , the suppression is a factor of 0.84. 

The separation between the maximal CP violating case 
8 = \ and the CP conserving case 8 = at the end of 
the fcth period is 



V —V 



1 



'13 sin 2023 (sin 7 + cos 7) sin 2kQ. 

(54) 

The maximum separation between these two curves oc- 
curs whenever £;$ = ^ + ~ and 7 = J , ~ ; this separa- 
tion is 



V. 



< 



1 

71' 



013 sin 20 



23 



(55) 



Unlike the separation between the 5 = 0, tt curves, the 
implementation of these constraints is nontrivial. For 




2L 3L 
Baseline (L) 

FIG. 5: (Color online) Oscillation probability v e — > 
through a castle wall potential using the same oscillation pa- 
rameters as in Fig.Hbut with E = 200 MeV, pA = 3.8 g/cm 3 , 
L A = 4024 km, p B = 11.5 g/cm 3 , L B = 924 km. The curves 
are as in Fig. [4] 



k = 1, the requirement for $ implies g = sin 2 $ = |X| 2 . 

With A3 = and the requirement on 7, these constraints 

^2 



translate into Xf = A = A|. For fixed vacuum values 
of the mixing angle 0i2 and mass-squared difference A21, 
there are five remaining free parameters for a general cas- 
tle wall profile: the neutrino energy and the density and 
length of the two regions in the castle wall. With the 
three constraints on the values of Xj , this leaves two free 
parameters, say, the two densities; however, for general 
values of the two densities, applying these constraints 
can result in a complex value for the neutrino energy. 
Depending on the value of 0i2, the densities pa and pb 
must be sufficiently different in order for the constraints 
to result in a real value for E. In Fig. [5] we show an 
example of these constraints which maximizes the sepa- 
ration between the 6 = [black] solid and 8 = % [green] 
dotted curves at the end of each odd period. The param- 
eters which we employ to produce these curves results in 
the phase values <I> = 37r/4 and 7 = 7r/4. If one knew the 
values of all the neutrino oscillation parameters, save the 
CP phase, this castle wall would could resolve 6 = and 
6 = I7, but a degeneracy in the parameter space would 
remain if measurements were made only at the bound 
aries because the 8 



0, ^ curves intersect as do the 



5 = f,7T curves. 

As an aside, we comment upon one additional case of 
parametric resonance which carries the phase condition 
7 = 0, 7r; as in the previous example, this extra condition 
places severe restrictions on the baseline and neutrino 
energies. From the definition of the phase 7 in Eq. (30 1, 
it is clear that A2 must vanish when 7 = 0, n. Thus, 
along with the condition for parametric resonance, one 
requires Ai = 1. To satisfy all of these requirements, in 
one region, say region A, the width of the region La must 
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be an integer number of wavelengths so that sa = 0. In 
the other region, the width of the region Lg must be an 
integer-plus-one-half wavelengths so that cb = 0, and the 
energy of the neutrinos must be at the MSW resonance 
in region B so that c-iq b = 0. In spirit, these criteria do 
not reflect true parametric resonance; rather this is essen- 
tially a manifestation of the MSW resonance. Neutrinos 
which travel an integer number of wavelengths through 
region A will exit the region in essentially the same state 
as they entered the region; then, in region B, we have ex- 
plicitly required the energy to be at the MSW resonant 
value. Turning to the conditions require $ = J, %- so 
that it is not possible to simultaneously satifsy 7 = 0, tt 
and fc$ = ^f- + j ; thus, the curves for all values of 5 will 
intersect at the boundary between periods. In Fig.[6j we 
implement this scenario. The neutrino energy is chosen 
to match the MSW resonance of region A, and La is 
one-half wavelength in size. For the next region, Lg is 
chosen to be a full wavelength in size. For the different 
values of the CP phase, we see that the oscillation prob- 
ability is the same value, but that the curves are quite 
distinguishable in the interior of region B. 



DISCUSSION 

For neutrino energies on the order of < 1 GeV, the 
wavelength of vacuum neutrino oscillations associated 
with the mass-squared difference A21 is 3x 10 4 km. Given 
this scale, a closed laboratory demonstration of paramet- 
ric resonance is not feasible. The only recourse is to use 
the Earth's mantle-core-mantle transition to induce the 
resonance. Unfortunately, this means that the neutrinos 
can only travel through fewer than two periods of a castle 
wall potential; however, sub-GeV atmospheric neutrinos 



FIG. 7: (Color online) We plot (a) X3, (b) $, and (c) 7 versus 
the cosine of the zenith angle for neutrinos traveling along a 
chord through the earth. For the solid [black] curves, the 
neutrino energy is 200 MeV; for the dashed [red] curves, 400 
MeV; for the dotted [blue] curves, 800 MeV. 

We model the earth as a constant density core of radius 
R c = 3485 km and density p c = 11.5 g/cm 3 surrounded 
by a constant density mantle with p m = 4.5 g/cm 3 . The 
path of a neutrino, a chord through the earth's interior, 
can be parametrized via the zenith angle 0. Neutrinos 
which traverse the entire diameter of the earth would 
have = 7r and thus cos = — 1. For this trajectory, the 
amount of the mantle seen by the neutrino before enter- 
ing the core is L m — i? e — R c — 2886 km where the radius 
of the earth is R e = 6371 km, and the neutrino traverses 
a path through the core of length L c = 2R C = 6970 km. 
As decreases, the mantle path length increases while 
the core path length decreases; generally, one has 

L m = -R e cos - yjR 2 c - {R e sin 0) 2 , (56) 
L c = 2^Rl-{R e sin©) 2 . (57) 

At the zenith angle cr it, the neutrino's trajectory is tan- 
gential to the core; this occurs when sin0 crit = R c /R e ; 
i.e., cr ;t = 2.56, or cos0 cr it = —0.84. 

From these path lengths, we can compute the relevant 
data for parametric resonance as a function of 0. In 
Fig. [7J we plot X 3 , and 7, respectively, for select neu- 
trino energies from 200 MeV to 800 MeV. From Fig.^a), 
we see that the condition for parametric resonance is sat- 
isfied in this energy range for chords with zenith angles 
such that —0.93 < cos© < —0.9; this is consistent with 
the neutrino oscillograms in Refs. [5D1 [ST]. The para- 
metric resonance is reflected in the enhancement of the 
v e — ► oscillation probability that would be seen at a 
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detector located at the end of the chord after the neu- 
trino has traveled a total distance of 2L m + L c . In Fig.|8j 
this enhancement is apparent as we plot, as a function 
of cos 9, the detector value of T efl for 200 MeV and 800 
MeV neutrinos for various values of the CP phase 6. The 
average location of the peak values of T ef i for the differ- 
ent values of 8 corresponds to the zenith angle at which 
X 3 = 0. Though the overall amplitude of the 800 MeV 
curves are suppressed because the energy is far from the 
MSW resonance, it is noteworthy that there is a large sep- 
aration between the peak values of the oscillation prob- 
ability for different values of 8. 




-0.95 -0.9 

COS0 



-0.85 



FIG. 8: (Color online) The oscillation probability "P efl for a 
neutrino which has traveled along a chord through the earth 
with zenith angle O for energies (a) 200 MeV and (b) 800 
MeV. The solid [black] curves have 5 = 0; the dashed [red] 



curves have 8 = 7r; the dotted [green] 
dot-dashed [blue] curves have 5 = t? 



curves have 8 



the 



To explore this point further, we plot in Figs. [9] and 



10 the v e — > v„ oscillation probability for the 200 MeV 
and 800 MeV neutrinos for paths along the chords which 
correspond to parametric resonance, cosO = —0.905 and 
—0.919, respectively. As the terminus of the neutrino's 
path through the earth is not an integer number of pe- 



riods, we cannot use Eq. (49) to ascertain V eil here. At 
best, the previous analysis only informs our knowledge of 
the state as the neutrino leaves the core at L = L m + L r \ 



from Eq. (46 1, the state, in the propagation basis, is 



v{L) 



cos4> 
— ie 11 sin$ 



(58) 



where the phases <E> and 7 can be read off of the plots 
111 Fig. [7] For the remaining bit of the path through 
the mantle, one can simply operate on this state wi th a 
constant density evolution operator U A (L m ), Eq. (42), 
to determine the state at the end of the chord. 

Though the curves in Figs. [9] and [10] are rather dif- 
ferent, they can be understood within the same frame- 
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FIG. 9: (Color online) Oscillation probability v e — s> for a 
chord through the earth with cos O — —0.905 which corre- 
sponds to Lm = 3576 km and L c = 4376 km. The various 
curves correspond to different values of the CP phase 8 as in 
Fig-H 
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FIG. 10: (Color online) Oscillation probability — \ for 
a chord through the earth with cosO = —0.919 which corre- 
sponds to Lm = 3438 km and L c = 4834 km. The various 
curves correspond to different values of the CP phase 8 as in 
Fig. [4] 



work as both systems approximately satisfy the half- 
wavelength condition with smip m = 1 = siny c - A conse- 
quence of the half- wavelength is the fact that the phase 7 
takes the value ±f . For both the 200 MeV and 800 MeV 
trajectories, this phase is actually 7 ~ — 0.47r; however, 
we will assume 7 = — 7r/2. Given this, we deduce from 
Eq. Q that sin$ = sin(26» c - 29 m ). With this simplifi- 
cation along with the assumption for 7, we find that for 
the half-wavelength condition the oscillation probability 
at the end of the neutrino's trajectory is 



V ell {2L m +L c ) w c 2 23 sin 2 (48 m -26 c ) 

40 c ). (59) 

As 800 MeV is roughly eight times the MSW resonance 
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energy in the mantle, the overall amplitude of these os- 
cillations is suppressed relative to the 200 MeV case. So 
far as the half-wavelength condition is satisfied, the CP 
phase only enters in the second term of the righthand 
side of the previous equation. As it enters only as cos 8, 
there will be no difference between curves with 5 = f , Tp. 
The difference between the CP conserving and maximally 
violating cases is given by 



TWo " « 013S23C23 sin(80 m - 40 c ). (60) 



For a particular matter profile, the size of this difference 
is controlled by the factor involving the effective #12 mix- 
ing angles in the mantle and core. For the 200 MeV case, 
one has sin(80 m - A9 C ) « 0.31, and for the 800 MeV neu- 
trinos, this is sin(80 m — 46 c ) sa —0.81. This accounts for 
the large separation between the curves in Fig. [TUJ rela- 
tive to the 200 MeV case. In terms of maximizing this 
difference, | sin(86» m - 40 c )| = 1 for E = 500 MeV, so in 
principle, this energy would be best to differentiate the 
CP conserving or violating cases for atmospheric neu- 
trinos traveling through the earth along a chord which 
satisfies the half-wavelength condition. 

In an actual experiment, the finite energy resolution of 
the detector must be considered; fortunately, the effects 
discussed herein are not washed out by a broad spec- 
trum neutrino source. Focusing upon the optimal en- 
ergy to ascertain CP violation, 500 MeV, we consider a 
neutrino beam with a flat energy spectrum between 400 
MeV and 600 MeV. In Fig. [TTJ we plot the beam's oscil- 
lation probability V elJt along a chord through the earth 
parametrized by the zenith angle 0. The curves which 
have conserved CP symmetry, 8 — 0, it, and maximal CP 
violation, 8 — § , , should be experimentally discernible 
for this broadband source. 

The bulk of current data on atmospheric neutrinos 
comes from the Super-Kamiokande collaboration [4"3Tl4"6"] . 
Unfortunately, this data cannot be used to search for the 
effects discussed in this paper because, most crucially, a 
water Cerenkov detector cannot distinguish between neu- 
trino and anti-neutrino events and, secondarily, the abil- 
ity to correctly ascertain the incident neutrino's zenith 
angle from the (detected) charged lepton is poor below 1 
GeV. A proposed detector, the magnetized iron calorime- 
ter (ICAL), at the India Neutrino Observatory (INO) 
[47) can distinguish muon neutrinos from anti-neutrinos, 
overcoming the primary impediment of a water Cerenkov 
detector. The resolution in zenith angle for the ICAL is 
still poor at low energies [48] . If we were to include this 
additional uncertainty in our computations, the peaks in 
Fig. [TT] will become smeared out, but the curves should 
still be experimentally distinguishable. 
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FIG. 11: (Color online) The oscillation probability V ell for a 
neutrino which has traveled along a chord through the earth 
with zenith angle O averaged over a flat neutrino spectrum 
from 400 MeV to 600 MeV. The various curves correspond to 
different values of the CP phase 8 as in Fig. [8] 

CONCLUSION 

Using an approximation appropriate for sub- GeV neu- 
trinos traveling through the earth, we are able to study 
in a three-neutrino framework parametric resonance of 
neutrino oscillations for a periodic density profile. Com- 
mensurate with the initial level of approximation, we de- 
velop a parametric resonance condition similar to the ex- 
act condition for two-neutrino systems. For a castle wall 
density profile, it is shown that the v e — > za, oscillation 
probability at an integer number of periods is enhanced 
and bounded by cos 2 #23. The CP phase 8 enters into 
the oscillation probability at these points via a phase ip 
which is proportional to #13 and involves the phase 7, a 
characteristic of the density profile and neutrino energy. 
This phase is present in an exact two-neutrino framework 
of parametric resonance but is of no measurable conse- 
quence in that context. 

As expected, in the three-neutrino framework paramet- 
ric resonance significantly enhances the oscillation prob- 
ability. We examine in detail instances of parametric 
resonance in which the phase 7 takes on three different 
values. When 7 = f , the v e — > oscillation proba- 
bilities achieve maximal separation at the end the first 
period for S = and 8 — n, provided <& = J. When 
7 = ? , this oscillation probability achieves maximal sep- 
aration for the 8 = and 8 = f cases, given the same 
condition for $. Though these trajectories are best for 
differentiating 8 = and 8 = f, they also suffer degen- 
eracies for 8 = 0, Tp and 8 = §,7r. Finally, whenever 
7 = and $ = |, the V eli oscillation probability for all 
values of 8 is the same an the end of each period. 

We also apply this formalism to sub-GeV neutrinos 
which travel along a chord through the earth, using the 
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mantle-core transition to generate parametric resonance. 
Significant enhancement of the oscillation probability ex- 
ists even in the case in which the neutrino energy is far 
from the MSW resonance. Though a path through the 
earth is not an integer number of periods, the formal- 
ism is useful to determine the state of the neutrino upon 
leaving the core. As the trajectories through the earth 
nearly satisfy the half-wavelength condition, the oscilla- 
tion formulae simplify greatly. Insofar as this condition 
is satisfied, we note that energies near 500 GcV will be 
best for differentiating the S = and 5 = f cases in the 
v e — > oscillation channel. 

As for a clean experimental confirmation of these res- 
onant oscillations, two main impediments exist. First, 
the ability to differentiate between neutrino and anti- 
neutrino events is crucial. Water Cerenkov detectors like 
Super-K do not have this ability, but ICAL at INO could 
overcome this. The second issue deals with the ability to 
accurately assess a neutrino's incoming zenith angle O. 
For sub-GeV neutrinos, the ICAL detector has poor res- 
olution for the zenith angle, but at 500 GeV, a detector 
should still be able to distinguish the cases of S = 0, ir and 
S = §, ■ Regardless, a semi-analytic understanding of 
the interplay between parametric resonances and CP vi- 
olation for sub-GeV neutrinos provides a useful backdrop 
for future data analysis. 
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